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Abstract 

We present, using the statistical model, a new analysis of hadron production in 
central collisions of heavy nuclei. This study is motivated by the availability of 
final measurements both for the SPS (beam energies 20-160 AGeV) and for the 
RHIC energies (y'ijv7V=130 and 200 GeV) and by updates in the hadron mass 
spectrum, which is a crucial input for statistical models. Extending previous studies 
by inclusion of very high-mass resonances (m> 2 GeV) , and the up-to-now neglected 
scalar a meson leads to an improved description of the data. In particular, the 
hitherto poorly reproduced energy dependence of the /t^^ ratio at SPS energies 
("the horn") is now well described through the connection to the hadronic mass 
spectrum and, implicitly, Hagedorn's limiting temperature. We thereby demonstrate 
the intimate connection between the horn and the QCD phase transition. 



1 Introduction 



One of the major goals of ultrarelativistic nuclear collision studies is to obtain information 
on the QCD phase diagram [1]. A promising approach is the investigation of hadron 
production. Hadron yields measured in central heavy ion collisions from AGS up to RHIC 
energies can be described very well P|5|HI5|S|7|I51I51[TI]|11|12II13|14] within a hadro-chemical 
equilibrium model. In our approach [2]|^1[71[T2] the only parameters are the chemical freeze- 
out temperature T and the baryo-chemical potential /ij, (and the fireball volume V , in case 
yields rather than ratios of yields are fitted); for a review see [T5] . 

The main result of these investigations was that the extracted temperature values rise 
rather sharply from low energies on towards -^/snn —10 GeV and reach afterwards con- 
stant values near T=160 MeV, while the baryochemical potential decreases smoothly as 
a function of energy. The limiting temperature [16] behavior suggests a connection to the 
phase boundary and it was, indeed, argued [17] that the quark-hadron phase transition 
drives the equilibration dynamically, at least for SPS energies and above. Considering 



also the results obtained for elementary collisions, where similar analyses of hadron mul- 
tiplicities, albeit with several additional, non-statistical parameters (see [T81IT9] and refs. 
therein), yield also temperature values in the range of 160 MeV, alternative interpreta- 
tions were put forward. These include conjectures that the thermodynamical state is not 
reached by dynamical equilibration among constituents but rather is a generic fingerprint 
of hadronization [20121] . or is a feature of the excited QCD vacuum [22]. The analysis 
presented below lends further support to the interpretation that the phase boundary is 
reflected in features of the hadron yields. 

While in general all hadron yields are described rather quantitatively, a notable exception 
was up-to-now the energy dependence of the ji:^ ratio which exhibits a rather marked 
maximum, "the horn", near y'J/vjv — 10 GeV [23]. The existence of such a maximum 
was, in fact, predicted [21] within the framework of the statistical model, but the ob- 
served rather sharp structure could not be reproduced [12]. Other attempts to describe 
the energy dependence of the j-n'^ ratio within the thermal model [TTfT^ were also 
not successful, except when an energy-dependent light quark fugacity was used as an ad- 
ditional parameter [T3]. Furthermore, all attempts to reproduce this structure within the 
framework of hadronic cascade models also failed, as is discussed in detail in [23]. As a 
consequence, the horn structure is taken in [23] as experimental evidence for the onset of 
deconfinement and quark-gluon plasma formation, and as support for the predictions of 

In this letter, we present a new analysis of hadron production in central nucleus-nucleus 
collisions. The motivation for this is twofold: i) on the experimental side the data set has 
been recently extended (and/or finalized), in particular for the SPS energy regime by the 
NA49 collaboration [23|26|27j and for the RHIC energies by the STAR collaboration [21] , 
and ii) we would like to explore the consequences of an improved hadronic mass spectrum 
in which the a meson and many higher-lying resonances are included. We note that, with 
the exception of the a meson, the full hadronic mass spectrum has already been used in 
our recent investigation [18] of hadron production in e^e~ collisions. In the following we 
first discuss the update in the hadronic mass spectrum and then explore its consequences 
for the description of all available data from SIS to RHIC energies. 

For the hadronic mass spectrum we are using the complete mass states published recently 
by the Particle Data Group (PDG) [29]. Not all the states have known branching ratios 
and in such cases we have assigned values based on analogies to the nearest states with the 
same quantum numbers. The systematic uncertainty introduced by this is estimated to be 
well below the bias introduced if the high-mass states would not be considered. Another 
recent development in the field of hadron spectroscopy is the strengthening of the case for 
the existence of the a meson (labelled /o(600) in the PDG compilation [29]) [30|31|32j . The 
(J meson is a broad structure, whose properties are extracted from fits of measurements 
in various channels (see [30] for a recent review) and decays into vr+vr^. We have adopted 
as "nominal" the values for the mass and (Breit-Wigner) widthQ as extracted in ref. [30] , 
mo-=484 MeV, rcr=510 MeV and will investigate the effect of different parameters, namely 
also the case m^=600 MeV, r<,=600 MeV. 



We ignore here that the width of the a meson is not small compared to its mass. Further 
investigations will have to deal with this issue but we note that for the p meson the situation is 
not so different. 
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Figure 1. The energy dependence of the increase in pion yield after inclusion of high- mass 
resonances (dashed) and the a meson, characterized by its mass and width (m„,T(j) in MeV. The 
dash-dotted line depicts the energy dependence of the temperature [12) used in the calculations. 

To explore the consequences of the improved mass spectrum we show as a function of 
energy, in Fig. [1], the increase of yields of pions after inclusion of high-mass resonances, 
relative to the case when hadrons up to a mass of 2 GeV were considered, as in our earlier 
study [12], and after the inclusion of the a meson. The observed energy dependence is 
driven mostly by the change in temperature, also shown in Fig. [H The parametrizations 
for the energy dependence of T and fib established in [12] were employed. The high-mass 
resonances lead to an increase of about 13% for the calculated pion yields. This increase 
levels off near the point where the temperature reaches its limiting value, thereby sharp- 
ening the structure in the K~^/7r~^ ratio, as will be demonstrated below. A further 3.5% 
increase of the calculated pion yield is due to the presence of the a meson, with a rather 
small dependence on its mass and width. For 771^=600 MeV, rCT=600 MeV the calculations 
lead to about 1% fewer pions. For the A hyperons, the new high mass resonances lead to 
an increase in the calculated production of about 22%, while the addition of the a meson 
has no effect. An increase of up to 6% is observed for protons, while for kaons this increase 
is about 7%. 

We have shown earlier [12] that the thermal fits of hadron yields and of ratios of yields 
lead to very similar results. For the present analysis we focus on fits of yields. We mostly 
utilize mid-rapidity data, but, at lower SPS energies, fit also the hadron yields for the full 
phase space. 

In Fig.[2]we present a comparison of measured and calculated hadron yields at the energies 
of ^/snn=7.6 GeV (beam energy of 30 AGeV at SPS) and y^J/v7v=200 GeV. The model 
is successful in reproducing the measurements and this applies to all energies, from 2 
AGeV beam energy (fixed target) up to the top RHIC energy of ^/sJ^=200 GeV. The 
reduced values are reasonable. In most cases the fit quality is improved compared to 
our earlier analysis |T2], even though the experimental errors are now smaller. Whenever 
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Figure 2. Experimental hadron yields and model calculations for the parameters of the best fit 
at the energies of 7.6 (left panel) and 200 GeV (right panel; the 0, yield includes both 0" and 

n+). 

several independent measurements are available, we have employed a weighted mean of 
the data following the recipe given in the introduction of [29]. A special case is that of 
the top SPS energy {y/sNN=17.3 GeV), where the disagreement between the NA49 and 
the NA57 data persists. As in the case of our earlier analysis [12], we have moved the 
difference in the fit results into the respective systematic error. A disagreement between 
the experiments is seen at the top RHIC energy for pions and protons, see Fig. [2l which 
is the reason of the large reduced x^- A fit of ratios is in this case more suited, but we 
note that a fit of the STAR yields alone gives T=162 MeV, /ife=32 MeV, V=2A00 fm^, 
with a very good x^/A'^[i/=9.0/ll. The resonances were not included in the fits, but are 
quite well reproduced by the model. 

An important result of our analysis is that the resulting thermal parameters are close to 
those obtained earlier [12] and are in agreement with other recent studies [T3|28] . This 
indeed confirms that the common practice of including in the thermal codes hadrons up 
to masses of 2 GeV (for instance in the publicly- available code THERMITS [33]) does 
not lead to significantly biased fit parameters. Nevertheless, there are small variations. In 
Fig. [3] we present the energy dependence of T and fib in comparison to our earlier results 
. We have parametrized T as a function of ^/sNN with the following expressioiH]: 



T, 



lim ' 



1 + exp(2.60 - ln(./sivw(GeV))/0.45) 



with the "limiting" temperature Tiim=^Q4: MeV. This value is slightly higher compared 
to our earlier value of 161±4 MeV [12] due to the higher temperatures presently derived 
for the RHIC energies. The approach to is presently more gradual compared to our 
earlier parametrization. 

The values of /if, extracted for the two lowest SPS energies deviate somewhat from the 
continuous trend suggested by all the other points. At these energies the fit to data in full 



For fib, there is no need to change our earlier |T2] parametrization: /if, [MeV] 
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Figure 3. The energy dependence of temperature and baryon chemical potential at chemical 
freeze-out. The results obtained here are compared to the values obtained in our earlier study 
|12j . The lines are parametrizations for T and Hf, (see text). 

phase space does lead, as expected, to larger values of fib, which do fit in the systematics. 
At 40 AGeV (y^J/v7v=8.8 GeV) the resulting values of T and /if, from the fit of An data 
are very similar to those obtained from midrapidity data. 

We employ the above parametrization to investigate the energy dependence of the relative 
production yields /tt^ and A/vr", shown in Fig. HI The K'^ jix'^ ratio shows a rather 
pronounced maximum at a beam energy of 30 AGeV [23], and the data are well repro- 
duced by the model calculations. In the thermal model this maximum occurs naturally at 
^ysNN —8 GeV |2l]. It is due to the counteracting effects of the steep rise and saturation 
of T and the strong monotonous decrease in /if,. The competing effects are most promi- 
nently reflected in the energy dependence of the A hyperon to pion ratio (lower panel of 
Fig. Hj), which shows a pronounced maximum at ^/sNN —5 GeV. This is reflected in the 
/n^ ratio somewhat less directly; it appears mainly as a consequence of strangeness 
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Figure 4. Energy dependence of the relative production ratios j-K^ and K/tt" . With the 
dotted hne we show for the /tt'^ ratio an estimate of the effect of higher mass resonances (see 
text). The dashed hnes show the energy dependence of T (upper panel) and /i;, (lower panel). 

neutralitjll], assumed in our calculations. 

The model describes the /n^ data very welH over the full energy range, as a con- 
sequence of the inclusion in the code of the high-mass resonances and of the a meson, 
while our earlier calculations [12] were overpredicting the SPS data. At RHIC energies, 
the quality of the present fits is essentially unchanged compared to [I2], as also the data 
have changed somewhat. The model also describes accurately the A/tt" data. We note 
that the maxima in the two production ratios are located at different energies [35]. The 
model calculations reproduce this feature in detail. 

The calculated /ti^ ratio (solid line in Fig. Hj) is likely to decrease further at energies 
beyond the maximum and the peak is likely to sharpen somewhat if our presumably 

^ Recent studies within the UrQMD model [M] suggest the possible presence of net strangeness 
at midrapidity; the present analysis lends no support to this conjecture. 
^ On a statistical basis, with a relatively good x^/-^d/=21/12. 



still incomplete knowledge of the hadronic spectrum for masses larger than 2 GeV would 
improve. To get a qualitative idea about the size of the effect we have assumed that the 
true mass spectrum rises exponentially, with a Hagedorn temperature parameter Th ~ 
200 MeV (as observed in the mass range up to 2 GeV). Based on this spectrum we have 
calculated the modification of the K"^ /n~^ ratio by inclusion and exclusion of states above 
3 GeV. The strongest contributions from these states to kaons comes from the decay 
of K* mesons (where only 1 kaon is produced for each additional K*) while all high 
mass resonances produce multiple pions, thereby reducing the /ti^ ratio. Using these 
assumptions and Boltzmann suppression of high masses we arrive at the dotted line in 
Fig. m While the qualitative nature of the estimate should be underlined, it shows in any 
case the expected features, namely a sharpening of the peak and a further reduction in 
at high energies. The uncertainty of the calculations due to the mass and width 
of the a meson are at the percent level only. Another few percent uncertainty, which is 
difficult to assess quantitatively, arises from the unknown branching ratios of the high- 
mass resonances. 



Table 1 

Predictions of the thermal model for selected hadron ratios in central Pb+Pb collisions at LHC. 
0.072 0.164 0.163 0.042 0.0054 0.00093 



In Table [T] we present our updated predictions for the LHC energy {^JsNN='^>■'^> TeV) [36] . 

In summary, we have demonstrated that by inclusion of the a meson and many higher mass 
resonances into the resonance spectrum employed in the statistical model calculations an 
improved description is obtained of hadron production in central nucleus- nucleus collisions 
at ultra-relativistic energies. The most dramatic improvement is visible for the /n^ 
ratio, which is now well described at all energies. The "horn" finds herewith a natural 
explanation which is, however, deeply rooted in and connected to detailed features of 
the hadronic mass spectrum which leads to a limiting temperature and contains the QCD 
phase transition [TB]. It is interesting to note that central questions in hadron spectroscopy 
such as the existence (and nature) of the a meson apparently play an important role in 
quark-gluon plasma physics. Our results strongly imply that hadronic observables near 
and above the horn structure at a beam energy of 30 AGeV provide a link to the QCD 
phase transition. Open questions are whether the chemical freeze-out curve below the horn 
energy actually traces the QCD phase boundary at large values of chemical potential or 
whether chemical freeze-out in this energy range is infiuenced by exotic new phases such 
as have been predicted in [37j. In any case these are exciting prospects for new physics 
to be explored at RHIC low energy runs and, in particular, at the high luminosity FAIR 
facility. It remains a challenge, though, to identify clear signals for quark-gluon plasma 
properties in this energy range. At the high energy frontier, the measurements at LHC 
will be a crucial test of the present picture. 
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